A Josephson junction as a detector of Poissonian charge injection 
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We propose a scheme of measuring the non- Gaussian character of noise by a hysteretic Josephson 
junction in the macroscopic quantum tunnelling (MQT) regime. We model the detector as an 
(under) damped LC resonator. It transforms Poissonian charge injection into current through the 
detector, which samples the injection statistics over a floating time window of length ~ Q/ooj, where 
Q is the quality factor of the resonator and ouj its resonance frequency. This scheme ought to reveal 
the Poisson character of charge injection in a detector with realistic parameters. 

PACS numbers: 72.70.+m,73.23.-b,05.40.-a 



Presently there is considerable effort on characteriz- 
ing and measuring the statistics of electrical current and 
the non-Gaussian nature of its fluctuations in mesoscopic 
conductors PJQ- The discrete nature of charge injection, 
e.g., in tunnel junctions, with typically Poissonian statis- 
tics can be revealed by studying not only the average 
current (I) and its variance (SI 2 ), unlike in the case of 
normally distributed current, but higher moments also, 
most notably the third (central) moment (SI 3 ). Nor- 
mally these higher moments introduce very small signals 
and the filtering requirements are strict, because of which 
it is very hard to measure them (see, e.g., |UH). There 
is, however, one recent experiment which succeeded to 
demonstrate the existence of non- vanishing third moment 
in transport through a non-superconducting tunnel junc- 
tion Q. The measurement required very long averaging 
times. Therefore, alternative techniques to collect more 
information, and perhaps eventually to determine further 
higher moments, are definitely needed. In their recent ar- 
ticle Tobiska and Nazarov p proposed an overdamped 
Josephson tunnel junction (array) as a threshold detector 
to measure such full counting statistics (FCS) making use 
of rare over-the-barrier jumps arising from current fluctu- 
ations. The nearby on-chip detector is a definite benefit 
of this proposal due to its natural high bandwidth. Yet 
they considered the limit where tunnelling is perfectly 
suppressed whereby the set-up becomes experimentally 
less accessible. In our proposal we consider an under- 
damped single Josephson junction (or a DC-SQUID) in 
the MQT limit. We demonstrate that the experimental 
complications of the proposal [5|, e.g., the need of a mul- 
tij unction Josephson junction array, and the fact that the 
escape threshold is more difficult to measure in an over- 
damped junction, are overcome in our scheme, and show 
that the effect of higher moments is pronounced using a 
threshold detector with parameters deduced from earlier 
MQT experiments (see, e.g., Ref. |6j). 

We discuss a simplified model where charges are in- 
jected according to Poisson statistics on a Josephson 
junction, which in turn is described by a damped har- 
monic oscillator [LCR, or a linearized resistively and ca- 
pacitively shunted junction (RCSJ) model]. We show 
how this environment performs a conversion from dis- 
crete (charge) statistics into continuous (current) statis- 




FIG. 1: The scheme of the threshold measurement and the 
model used to extract the results of this paper. The current 
through the Josephson junction, depicted as a DC-SQUID, is 
determined by the currents through the scatterer (I±) and the 
ideal bias line (h), respectively. Switching from the super- 
current state to the normal state of the Josephson junction is 
signalled by a non-zero voltage V at the output of the ampli- 
fier. The resonator model of the Josephson junction and the 
charge injection are discussed in the text. 



tics. The proposed measuring scheme is sketched in Fig. 
There are two injecting lines, with currents I\ and 
I2. I\ is generated by a voltage bias across the scatterer 
(for example a tunnel junction), and I2 runs in a directly 
connected line to be discussed below. We neglect the in- 
fluence of the transmission line connecting the injecting 
lines and the measuring Josephson junction. Capacitance 
of the injecting junction can be included in the RCSJ 
capacitance, and similarly the additional injection line 
(I2) and the connection to the voltage amplifier can be 
modelled as a parallel inductance, which may reduce the 
resolution of the detector as will be discussed at the end. 
For practical implementations we depict the detector as 
a DC-SQUID with critical current tunable by magnetic 
flux <£. 

The detector is driven by elementary charges e* 
"dropped" on the LCR resonator according to Pois- 
son statistics. Therefore, each event poses an elemen- 
tary current e*5(t — tj) into the resonator such that 
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equation whose solution it is: 
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FIG. 2: Simulated I(t) (thin solid lines) with the following 
parameters. / = 100, p = 0.01 and with two values of quality 
factor, Q = 1 and Q = 10, respectively. Thick solid lines are 
expectations of (/(£)) according to Eq. 



= ^Zj^5(t - tj) as illustrated in Fig. [I] Here 
tj is the time instant of the jth event and S(t) is the 
Dirac delta function. (Alternatively we could consider a 
voltage step of height e*/C at each instant tj.) Our task 
is to evaluate the current / and its moments through the 
junction (inductor L). The elementary current i(t,tj) 
through the Josephson junction at the time instant t, 
due to this jth event can be easily solved, and is given 
in the case of an underdamped resonator (quality factor 
Q = R/^L/C> 1/2) by 

i(t, tj) =e*^j(l - ^)~ 1/2 exp[-±uj(t - tj)] 

M^-^) 1,2 ^{t-tj))e{t-tj). 



Here, coj = (LC)~ X I 2 is the plasma frequency of the junc- 
tion and 0(t) is the Heaviside step function. Analogous 
results for overdamped case (Q < 1/2) exist, but they 
will not be considered explicitly here, since the detector 
is assumed to switch from supercurrent state to resistive 
state. Therefore we assume Q > 1/2 in what follows un- 
less otherwise noted. It is interesting to see some proper- 
ties of these elementary oscillations. At the time instant 
tj all charge e* is stored in the capacitor C, whereafter 
it performs damped oscillations at (angular) frequency 
ujj. On integrating ij(t) one finds a consistent result: 

h(t) induces current I(t) through the Josephson junc- 
tion that can be expressed as the sum of the elementary 
currents of Eq. on the basis of the linearity of the 



where now the time instants tj are perfectly uncorrelated. 
One can numerically simulate the current following Pois- 
son principle: divide time into very small intervals St and 
set the probability p <C 1 for one charge to tunnel during 
this interval. No multiple events in any of the St intervals 
are allowed. The average current 7, p and St are related 
through 7 = pe* /St. Figure [21 illustrates results of sim- 
ulation with two different values of Q, time unit is ooj 1 , 
and 7 = 100, in units e*cuj. The initial rise and oscil- 
lations are due to the fact that the injection of charges 
is suddenly initiated at t = 0. At larger values of Q the 
junction tends to oscillate at (angular) frequency cjj, but 
it is driven by random events and thereby dephased. 

Simulation yields a representative example based on 
the random numbers determining the time instants tj . It 
is more general and interesting to investigate the proper- 
ties of the distribution of current. Let us first consider the 
various moments of current I(t). This can be done in the 
following way. We take a long enough time interval r such 
that the interesting instant t at which we want to evaluate 
the moments of current satisfies t < r. Now we consider 
different ensembles of instants tj, such that N charges 
are injected within r, and then weight all these config- 
urations by the Poisson probability Pp isson(A^)- Since 
all the j = 1,...,7V events are uncorrelated and evenly 
distributed over < tj < r, we may write for the nth 
moment of I(t): 
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It is straightforward to integrate for (I n (t)) using i(t,tj) 
of Eq. QJ. Below we summarise results for the three 
lowest moments of I(t). The average current reads 

</(*))= 7{l-exp(-^) 
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Here we have identified / = e*(N)/r = 
— J^jvLi ^Ppoisson(A0, which is the asymptotic value of 
(I(t)) on t — » oo. The thick lines in Fig. [3 show the re- 
sult of (I(t)) using Eq. (@J), which indeed seem to follow 
the mean of the simulated curves as expected. In what 
follows, we drop out the argument t, and consider only 
results after the initial transient, i.e., t ^> 2Q/coj. The 
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second raw moment reads (I 2 ) = ^^p-I+ ^- T (N(N-1)), 
where (N(N - 1)) = En=i N ( N " l)ftoLon(W). The 
more interesting second central moment, the variance 
{SI 2 ) = ((J - (J)) 2 ), then reduces to 
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(5) 



After a straightforward derivation we similarly obtain the 
third central moment, (SI 3 ) = ((/ — (i)) 3 ), reading 



(si 3 ) = 



3(l + 2/Q 2 ) 



(e*u>j) 2 I. 



(6) 



According to Eq. 0, the shot noise of the injecting 
junction, 2e*7, is amplified by the quality factor of the 
resonator over the band whose width is ~ cjj, the (max- 
imum) response frequency of the detector. A measure of 
the non- Gaussian character of the current distribution is 
its skewness 7J, defined as S = (SI 3 ) / (5 1 2 ) 3 / 2 , which, 
according to Eqs. (|5J) and ® reads 



S = 



2 5/2 



3(Q 3 / 2 + 2Q-V2) 



(,3*^)1/27-1/2. 



(7) 



Results (|SJ)-([7| hold also in the overdamped case (Q < 
0.5). It is interesting to note some general features of 
S in Eq. J7J). The non-Gaussian "strength" increases, 
in accordance with the central limit theorem, with de- 
creasing I (less events recorded). The detector exhibits 
a memory of events over a time ~ 2Q/coj in the under- 
damped [l/(Qcoj) in the overdamped] case, and there- 
fore, the skewness attains its maximum value close to 
the crossover between underdamped and overdamped be- 
haviour, Q — 1: here the memory of the detector is short- 
est, and it responds to only a small number n of Poisson 
distributed events through the scatterer. In the example 
discussed below n ~ (I /e*)(2Q/uj) ~ 40. 

Figure [3| shows an example of the simulated (3-10 5 rep- 
etitions) current distribution at the time instant cojt = 40 
(far enough after the initial transient, cojt ^> 2Q) for 
the injected current with / = 10e*o;j and for a detector 
whose Q = 2. The main frame with logarithmic verti- 
cal scale shows the number of counts (in solid circles), 
demonstrating how the simulation differs from the Gaus- 
sian fit shown by the solid line. There are more hits 
at high currents in the simulation as compared to the 
normally distributed events as one would expect. The 
same is shown as the true distribution (counts divided 
by 3 • 10 5 ) on linear scale in the inset. Table [I] gives a 
comparison between the simulated values and the theo- 
retical predictions [Eqs. Q-©] of the three lowest mo- 
ments. The correspondence is satisfactory, although the 
variance (SI 2 ) falls outside the la uncertainty margin. 

Next we make a judgment of whether such a threshold 
detector provides a viable means to measure FCS. To 
this end we need to consider the escape rates from the 
super current state. We assume low temperature T such 
that thermally activated switching is suppressed. This 
is the case when T < hw j / '(27t/cb) , which is an easily 
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FIG. 3: Simulated current distribution at cujt — 40, Q = 2 
and / = 10e*o;j. The main frame shows the number of counts 
(out of 300 000) in solid dots, and the solid line is the best 
Gaussian fit to it. The inset shows the same results but as 
the true distribution, = counts/300 000, and on linear scale. 



TABLE I: Results 



simulation theory 

(/) 10.01=b 0.01 10.00 

(51 2 ) 9.91± 0.03 10.0 

(51 3 ) 4.65± 0.22 4.44 



accessible regime experimentally |g. Escape rate in the 
MQT regime in the presence of Gaussian noise has been 
discussed, e.g., in Ref. 9]. Here, we will allow for more 
general current statistics. Using the standard decay law 
we find that the probability of escape to the resistive 
state is given by 



P = 1 — exp 



[-/ 

J to 



to+At 



T(I(t))dt], 



(8) 



where At is the duration of the current pulse starting at 
t = to over which we monitor escape statistics. T is the 
current dependent escape rate in the MQT process for 
which one can find explicit expressions that depend on 
the junction and circuit parameters 0: r = iiexp(-B), 
where A = xV^J^U /(2ttH) and B = sAU/(hu;j). AU 
is the / dependent barrier height, and parameters \ an d 
s are Q-dependent. They assume values \ — 12^/^tt and 
s = 36/5 for large Q. For a pulse with At ^> Q/ooj 
we may write f^+ At T(I(t))dt ~ (r)At, where (T) = 

J-™ Here p(I) is the current distribution 

approximated, e.g., in Fig. 03 In the scheme of Fig. 
[^the average current I through the Josephson junction 
can be generated by any combination of the two (aver- 
age) currents I\ and I2 with the constraint I = T\ + 12. 
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FIG. 4: Escape probability P of a Josephson junction as a 
function of the average current I through the detector under 
different noise conditions. The shift between the solid and the 
dashed lines arises from Poisson statistics of charge injection. 
The dash-dotted line is the escape probability with the cor- 
responding Gaussian noise (same (5I 2 )). Escape in the ideal 
case of noiseless current is shown by the dotted line. 



Of particular interest are the cases where 1\ is either = 0, 
or it has positive or negative values of equal magnitude. 
Difference in the escape characteristics between the latter 
two cases provides a measure of the asymmetry of p(I) 
around its mean, the central topic of this article. Fig- 
ure 0] shows the escape histograms calculated under such 
three conditions using trapezoidal current pulses of du- 
ration At = 100 (is UUM We 

assume that the detector 
junction has a critical current Iq = 1 /aA [L = 7i/(2e7c)], 
C = 0.1 pF, and other parameters and p(I) are as in 
Fig. 03 (With these parameters At ^> Q/ujj.) We as- 
sume pure MQT escape with low dissipation [y,[l(j- The 
histograms are plotted as a function of average current 
7 through the detector driven by the two injection cur- 
rents in different proportions. The histogram shown in 



dotted line corresponds to no current fluctuations, i.e., 
all current is driven through the ideal line (7i = 0). The 
solid line is for the case when current through the scat- 
terer and that through the detector point in the same 
direction (|7| = |7i| + |^2|). With our circuit parame- 
ters I\ = +0.2/c- The dashed line is for the case when 
current 1\ points opposite to that through the detector 
(|7| = ||7i| - |7 2 || and 7i = -0.27c). The average shift of 
the 7i = ±0.2ic histograms with respect to the dotted 
line is due to the variance of current, whereas the pro- 
nounced shift between the last two is the more interesting 
effect of non-Gaussian current statistics. 

We conclude with a few practical remarks. The scheme 
presented here is simplified in many ways. Firstly, we 
do not take into account the (weak) dependence of the 
plasma frequency (Josephson inductance) on I to keep 
the discussion more phenomenological and transparent. 
This is fairly well justified for example in case of Fig. 21 
because all the relevant escape currents are smaller than 
0.8/co an d oc (1 — I / Iq) 1 ^ - Secondly, the escape his- 
tograms calculated assume low dissipation, which is not 
truly the case. The influence of dissipation on the MQT 
rate through environmental noise can be taken into ac- 
count by a minor scaling of x an d s parameters 9], and 
again to keep analysis on the basic level, we omit this 
since the effect is weak even when Q = 2. Thirdly, we 
assume that the injected charges do not produce current 
pulses either in the I2 current line, nor in the line to the 
voltage amplifier. This can be realised by large induc- 
tance (Lext) i n these lines, which in practice means long 
and narrow wires. If L ext ^> L, our argument is justified. 
Finally, the presented model is based on classical descrip- 
tion of the circuit dynamics: this is a valid starting point 
in the case of a tunnel junction scatterer whose tunnel 
resistance Rt > h/e 2 . 
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